Self-consistent calculation of electric potentials in Hall devices. 
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Using a first-principles classical many-body simulation of a Hall bar, we study the necessary 
conditions for the formation of the Hall potential: (i) Ohmic contacts with metallic reservoirs, (ii) 
electron-electron interactions, and (iii) confinement to a finite system. By propagating thousands of 
interacting electrons over million time-steps we capture the build-up of the self-consistent potential. 
The microscopic model sheds light on the the current injection process and directly links the Hall 
effect to specific boundary conditions at the particle reservoirs. 
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I. INTRODUCTION 



The calculation of self-consistent potentials of an in- 
teracting many-body system is of paramount importance 
for understanding the transport phenomena occurring in 
nanostructures and to develop realistic semiconductor 
device-models reaching from the classical to the quan- 
tum regime.— Ground state density functional theory 
is not always suitable to treat transport phenomena in 
semiconductors^ and time-dependent density functional 
theory is still limited to rather small systems (less than 
hundreds of electrons). On the other hand, classical the- 
ories for interacting electronic systems are well estab- 
lished, i.e. Thomas-Fermi screening and its extensions 3 * 
have been successfully used to obtain the effective mean- 
field potential of electrons in a semiconductor device. 
The presence of a stationary current through a device 
requires implementing an open-system, where the con- 
tacts provide sources and sinks for electrons. The con- 
tacts form an indispensable part of the simulated system. 
The methods described above are only partially able to 
address this situation, even for a purely classical system.— 
An interesting example is the classical Hall effect which 
gives rise to an unusual potential-theory problem^ fn a 
Hall bar, the voltages at the source and drain contacts 
are externally given and the current flow has to be cal- 
culated under the condition that the electrons do not 
cross the sample edges. These constraints on the current 
in the device and the voltages at the contacts have to 
be sa tisfied simultaneously. The resulting Hall potential 
(Fig. 1(b) I has been obtained via conformal mapping.^ 
The conformal map yields a unidirectional current flow, 
without counterpropagating currents along the edges. In- 
stead a diagonal current flow from one corner of the de- 
vice to the opposite corner emerges. Recently within 
the non-equillibrium network approach similar potentials 



were obtained 
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Mainly due to the lack of computational feasibility, a 
first-principles, microscopic model of the Hall effect has 
been lacking. But only such a model allows us to study 



and change important parameters (density, donor-layer 
distance, device geometry, gates, etc.) and to finally un- 
derstand which characteristics of a device, microscopic 
forces, and boundary conditions are necessary to gen- 
erate the experimentally observed Hall potentials^ To 
achieve this step, we have developed a new computational 
approach utilizing graphics processing units and adapted 
many-body algorithms used to study galaxy formation to 
nanodevices. In addition, our particle-based simulations 
highlight the dynamical nature of electron transport in a 
device, where we study the time-dependent build up of 
the self-consistent particle distribution. 

The manuscript is organized as follows: in Sect. HIl we 
describe the general idea of the particle-based simulation 
scheme and give the details of our method in Scct lllll Re- 
sults are discussed in Sects. |TV] |V]and finally we present 
conclusion and an outlook for future research in Sect. I VII 



II. 



INTERACTING PARTICLES SIMULATION 
OF SEMICONDUCTOR DEVICES 



Our first-principles model works in the spirit of molec- 
ular dynamics 1 - 1 - by using microscopic equations for the 
motion of electrons in a charge-field and highlights 
the importance of electron-electron interactions together 
with the boundary conditions at particle reservoirs. In 
the classical microscopic model, the force F k on the 
kth. electron consists of the forces due to all other 
N = N e + Nd electrons and donors, the confining poten- 
tial V w , and the velocity dependent Lorentz force F, 

F k = F^+Fl-VV w (r)\. 
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Ffc = q k r k x B, 

where B denotes the magnetic field. Each particle state- 
vector contains the position r k , velocity r^, and charge 
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<7fc (negative for electrons, positive for the stationary 
donors). In the following simulations we set the dielectric 
constant e = 8 and use an effective mass to* = 0.067 m e 
approximately matching electrons in a GaAs heterostruc- 
ture. 

Microscopic models of open systems present a recent 
development in theory^ Only the computationally very 
demanding particle-based simulations are able to accu- 
rately treat electronic transport on the classical level and 
to go beyond simplified drift-diffusion or hydrodynamic 
approaches^ Our simulations of large systems (10 3 — 10 4 
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FIG. 2: Left panel (color online): potential V(x,y,0) in the 
electronic subsystem S at z — due to the donor layer located 
at z = — 10 nm. Right panel: cut through the donor potential 
in the middle of the bar (see bright line in the left panel at 
x = 1.25 fim) showing oscillations with amplitudes of 2 meV 
due to the regular donor lattice. 
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FIG. 1: (color online) (a) System setup. The metallic reser- 
voirs R supply and remove electrons through the contact 
stripes C to/from the two-dimensional electronic subsystem 
S, which sits above the positively charged donor layer D. (b) 
Hall potential obtained from the conformal mapping method, 
where electrons enter the Hall device from the left source con- 
tact and drift to the right one. The contacts are metallic 
cquipotential surfaces. The black lines denote equipotential 
lines at —4,-2,0,2,4 mV. The dots at the left and right sides 
denote the positions of 500 contact points G used in the 
numerical simulation to model the Ohmic contact connect- 
ing the metallic reservoir with the electronic subsystem-layer, 
(c) Time-averaged potential obtained from the microscopic 
particle-simulation described in Sect. IIVI 



interacting classical particles) rely on a 200-1000 fold 
improvement of computational speed due to the use of 
high-performance graphics processing units, which allow 
us to simulate micrometer-size devices at realistic elec- 
tron densities. A quantum-mechanical calculation for 
such big open systems is still an impossible task. How- 
ever, classical simulations have always provided guid- 
ance for future quantum-mechanical simulations and are 
needed for understanding the quantum-to-classical tran- 
sition happening at higher temperatures. Especially 
for the understanding of the relation between classical 
and integer quantum Hall effect (IQHE) such simula- 
tions arc required, since often for simplified quantum- 
mechanical calculations the Fermi-liquid approximation 
is used, where interactions are effectively not present be- 
tween the quasi-particles. Notable exceptions are theo- 
ries of the fractional quantum Hall effect (FQHE) and 
the analysis of small systems with Coulomb-blockade 
physics^ Interestingly, electron-electron interactions are 
the corner-stone of the theory of the classical Hall effect 
and have to be considered in realistic semiconductor sim- 
ulations at higher temperatures X The particular shape of 
the classical Hall potential (Fig. 1) has also been observed 
in the IQHE (see Fig. 2 in RefEfl). 

Also in graphene devices the absence of the integer 
and fractional QHE in a four-terminal measurement in 
small samples is attributed to the influence of the Ohmic 
contacts, demonstrating the need to include the finite 
geometry, the contacts, and the interactions in the device 
model The theory of the IQHE in a finite graphene Hall 
device is discussed in RefUSI. Additionally, the IQHE 
displays a strong influence of the filling factor on the 
Hall potential, which is absent in classical models. 



III. BOUNDARY CONDITIONS, CONTACTS, 
AND RESERVOIRS IN OPEN SYSTEMS 



Our system Fig. |I(a) consists of several connected 
parts: two metallic three-dimensional electron reservoirs 
R(sourcc) and R(drain), two two-dimensional contact 
stripes C (source) over the complete device width between 
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x = [0,0.1] fim, and C(drain) between x = [2.4,2.5] (im 
where electrons from the reservoir are injected into a two- 
dimensional electron subsystem S. Uniformly distributed 
N c i = 8094 positive background charges are located in 
the planar donor layer D 10 nm beneath the electronic 
layer, whose extension exactly matches the extension of 
the electronic subsystem S. The resulting positive po- 
tential in the unfilled electronic layer (Fig. [2]) is bulged 
outward due to the finite extension of the system and 
has a range of 2.2-3.0 V. The motion of the electrons 
in the 2D subsystem is confined by a rectangular box- 
potential at y = 0, y = 1 /im, x = 0, and x = 2.5 /xm, 
which represents the etched borders of the nanodevice. 
Perpendicular to the plane, a homogeneous magnetic 
field B is present. The electrons sitting in the three- 
dimensional reservoirs are not part of our simulation and 
do not contribute to the potential calculation in the two- 
dimensional contact region, see Refill, model M4. The 
contact points Ci shown in Fig. 1(b) mark the positions 
of metallic spikes where electrons get transferred from 
the reservoirs to the subsystem S and vice versa. 

To find the self-consistent stationary electron flow 
through the system and the corresponding potential, we 
put initially 7800 electrons at random positions in the 
electronic system S. The initial distribution of the elec- 
trons was chosen to counterbalance the positive donor 
charges, but we do not enforce equal numbers of elec- 
trons and donors in the system, since we want to obtain 
the self-consistent solution. The reservoirs are connected 
to the contact points and enforce an equipotential sur- 
face within the contact stripes C. In a typical Hall mea- 
surement the source and drain reservoirs and contacts 
are kept on fixed potentials and one probes the Hall- 
voltage in the transverse direction along the y-axis. In 
mesoscopic physics, interactions are often considered to 
be absent in the semi-infinite leads connecting the de- 
vice region with the reservoirs. In our setup, we find 
that the electron-electron interactions within the con- 
tact stripes, from the moment of injection on, are cru- 
cial for the build-up of the self-consistent Hall-potential. 
We allow the number of electrons within the subsys- 
tem to change due to injection and removal events in 
the contact stripes C. Several techniques have been pro- 
posed to establish a microscopic model of carrier injec- 
tion from an Ohmic contact^^ Our introduction of con- 
tact points is motivated by studies of the microscopic 
nature of Ohmic contacts^ Ohmic contacts to quasi 
two-dimensional electron gases (2DEG) in AlGaAs/GaAs 
heterostructures consist of metallic spikes which pene- 
trate the cap-layer and reach down to the layer of the 
2DEG J£ In one dimension, the most realistic results have 
been obtained by using the reservoir method, where elec- 
trons can enter the contact stripe from a reservoir in or- 
der to maintain the applied bias-voltage locally in the 
contact stripe.— We have generalized the ID model to 
two dimensions and sample the potential in both con- 
tact stripes C at 500 contact-points Ci. We specify the 
target value for the respective contact potentials (here 
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FIG. 3: Time-evolution of the extremal value of the potential 
of all contact points in the source contact towards the target 
value V^om-cc = — 5 mV. The index i in this graph does change 
with time and refers always to the location of the contact 
point r(Ci,sourco) with the greatest deviation from the target 
value. 



Vsourcc = -5 mV, Vd ra i n = +5 mV). Initially, the elec- 
tron and donor charges produce an on (spatial) average 
flat potential over the device region, which does in gen- 
eral not match the target potentials at the contacts. 

Next we integrate the equations of motion of the in- 
teracting iV-body system in a magnetic field of B = 4 T 
with a time-step of At = 5 x 10~ 17 s using a modified 
Euler algorithm 
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r k (t + At) =r k (t) + At r k (t + At), 



where lui = -f^- denotes the Larmor frequency. Dur- 
ing the propagation we keep track of the time-averaged 
potential at all contact points Ci 
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where r(Ci) denotes the position of the contact point 
Ci. The starting point of the integration is initially set 
to to = 0. The integration loop gets interrupted each 
200At and we determine for source and drain separately 
the contact point with the maximum value of the time- 
averaged potential t^(C™o* rco ), V(Cf^ ain ) and the point 
with the minimum value V(C^J, V(C^J. How- 
ever, we discard the selected points in case the integra- 
tion time is smaller than a cyclotron period AT < T c in 
order to average out momentary potential fluctuations. 
Next we inject an electron at location r(C™o* rco ) in case 
t^(C™Q^ rcc ) > T4 ourcc and we mark the electron closest 
to r(Q™ rcc ) for removal if V(C^ mcc ) < V sourcc . The 
same process takes place at the drain contact. The rate 
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of 4 possible removals and additions every 200Ai corre- 
sponds to a maximum possible current of 16 fiA. After 
an injection or removal event, the lower limit to of the 
time-integration in Eq. [3J is set to the time of event. In 
our simulations we obtain currents of the order of 1 /iA, 
demonstrating that convergence was achieved and that 
there is no need to inject or remove electrons every 200At 
in order to maintain equipotcntial surfaces in the contact 
regions. Between injection and removal, the electrons 
move under the forces given by Eq. (JXJ) . 

The specific advantage of the used NVIDIA Tesla GPU 
board lies in the simultaneous calculation of 240 two- 
body interaction forces, which yields performances in the 
tcraFLOPS (10 12 floating point operations per second) 
ranged The limitation to single-precision arithmetic is 
of no major concern for the highly-chaotic classical sim- 
ulation, since we are interested in statistical averages. 
We perform 200 iterations steps on the GPU, before we 
transfer the instantaneous positions and velocities of all 
electrons to the CPU host. On the CPU, we calculate 
the potentials in the contact regions and inject and re- 
move particles according to the rules given above. Sub- 
sequently the updated set of particles is sent back to the 
GPU for propagating another 200 steps. The technical 
details of our hybrid GPU-CPU scheme will be described 
elsewhere. 

Finally, we calculate the resulting (time-averaged) po- 
tential at all points of the electron subsystem S. Our in- 
jection and removal procedure ensures that the potentials 
at the source and drain contacts converge towards their 
respective target values (see Fig. [3]), which are physically 
enforced by the presence of the metallic particle reser- 
voirs. If instead electrons are injected with equal prob- 
ability across the contact stripe, no Hall voltage across 
the device emerges. 
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FIG. 4: (color online) Cut through the time-averaged poten- 
tial across the center of the Hall bar at x = 1.25 fim. The 
target potentials are +5 mV (left source contact) and —5 mV 
(right drain contact). The oscillations in the potential profile 
are caused by the oscillations of the donor potential, shown in 
the right panel of Fig. [2l the blue circles denote the conformal 
map solution . Th e potential map of the complete Hall bar is 
given in Fig. 1(c) 



and the time- averaged value (see Fig. [4j). Figure [5] shows 
the histogram of the injection and removal events in 
the contact stripes, where for each bin we display only 
the net-result (number of injection events minus removal 
events) out of 10 5 events. The source contact predom- 
inantly injects electrons in the lower left corner, while 
the drain contact removes most electrons in the upper 
right corner. Thus the mean electron flow follows a di- 
agonal path across the device as displayed in Fig. |6(b)| 



IV. TIME-AVERAGED POTENTIALS AND 
ROLE OF INTERACTIONS 

We track the total number and the position of the elec- 
trons in the Hall bar. After a transient behavior (lasting 
ca. 0.1 ns) at the beginning of the simulation, a steady- 
state situation is reached where the electron number fluc- 
tuates around the value N e = 7975 ± 5, giving rise to an 
average electron density of n e = 3.19 x 10 15 m -2 . The 
current is obtained by introducing a flux-surface across 
the subsystem S or alternatively by counting the injec- 
tion and removal events at each contact i The net-current 
inflow from the source contact is +1.07 /iA, while the 
drain contact supports a net outflow of —1.07 /iA. The 
expected Hall voltage V xy 
related to n e via 



R xy I across the device is 



B 



7826 fi. 



(4) 



The calculated Hall voltage V xy = 8.4 mV is in good 
agreement with the target potential difference of 10 mV 
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FIG. 5: Spatial distribution of the net injection (positive axis) 
and removal (negative axis) events after a total of 10 events 
across the source and drain contacts. Electrons are predomi- 
nantly injected at the lower left corner of the source contact, 
and are most often removed at the upper right corner of the 
drain contact. 



The trajectory of an individual electron (Fig. |6(a) I can 
deviate considerably from the average behavior. Fig- 
ure [7] shows the distribution of the flight times of the 
electrons between the contacts. From the time-of-flight 
distribution we infer that most electrons are short lived 
and get absorbed within the same contact in which they 
got injected. These electrons do not contribute to the 
total current, but they are important to establish the 
mean Hall potential and the average current-density field. 
The average time-of-flight from source to drain contact is 
ttoi = 0.6 ns, and sets the time-scale required in order to 
obtain the converged result shown in Fig. 1(c) Starting 
the time-average at a later point t > i to { does not alter 
the picture. 

We also studied the importance of the electron-electron 
interactions by changing the dielectric constant e in a 
range from 1 — 8 and by additionally introducing a short- 
range cutoff of the Coulomb interactions. We find that 
for decreased Coulomb interactions (e > 10) the resulting 
potential does not resemble the conformal map result, 
but instead electrons start to pile up next to the con- 
verged contact region, while at the other side the electron 
density is reduced and results in a too positive value of 
the potential. Our finding shows that the electron must 
have a minimum degree of incompressibility in order to 
yield the classical Hall effect. Antisymmetrization and 
the Pauli principle (which are not part of the classical 
simulation) can provide other mechanisms to keep elec- 
trons apart and give an effective "Pauli incompressibil- 
ity". 



V. CONNECTION BETWEEN CLASSICAL 
AND QUANTUM HALL EFFECT 

As exemplified by our particle-simulation of a Hall- 
bar, the incorporation of the metallic boundary condi- 
tions together with the electron-electron interactions in 
the whole device (including the contact stripes) and the 
finite size of the system are all mandatory ingredients for 
the formation of a classical Hall potential. These facts 
were already noted by Hall in 1879^ p. 287, but the 
precise way how these ingredients lead to the Hall ef- 
fect could not been elucidated before due to the lack of 
computational methods. 

An intriguing problem is the passage from the clas- 
sical to the IQHE and the role of interactions in the 
IQHE. Recently it has been suggested that interactions 
are actually an important ingredient for the IQHE in 
order to calculate the critical exponents of percolation 
theories of the IQHE The experimental observation of 
the "classical" Hall potentiate under conditions where 
the resistivity is quantized provides another hint that 
electron interactions and the boundary conditions at the 
reservoir/subsystem interface plays an important role in 
the IQHE. Interactions enter the picture on two different 
length scales: while the Coulomb repulsion, the specific 
boundary conditions set by the device geometry, and the 
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FIG. 6: (a) Typical source-drain trajectory starting at the 
left source contact and reaching the right drain contact, (b) 
Time-averaged current density distribution. The arrow length 
is proportional to the local magnitude of the current density. 



metallic contacts are seen to enforce the global shape 
of the Hall potential, small clusters of neighboring elec- 
trons perform a highly correlated dancing pattern, which 
could be seen as the classical analogue to the FQHE3 
The interactions between the donors and the electrons re- 
sult in the screening of the positive donor charges, which 
become visible after performing the time-average of the 
instantaneous potential snapshots. 

A direct translation of the classical calculations to 
the quantum regime is not easily possible, since in 
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FIG. 7: Histogram of the time-of-flights of the injected elec- 
trons: number of source-to-drain paths (solid line), drain- 
to-source paths (dashed line), source-to-source paths (dotted 
line), and drain-to-drain paths (dash-dotted line). For com- 
parison, the fastest mean-field drift path yields a time-of-flight 
of 0.7 ns. 
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the quantum-case at low temperatures scattering events 
strongly depend on the occupation of initial states and 
the availability of final ones. Furthermore exchange and 
correlation effects have to be taken into account^! and 
the Coulomb interaction is modified due to the orbital 
extension of the electronic Landau levels. As a first 
step towards the incorporation of quantum-effects, we 
have repeated the simulations using a classical Hartree- 
approach, where Gaussian-shaped charge-clouds repre- 
sent the density of an electron within the first Landau 
level. The Gaussian distribution cuts off the Coulomb- 
interactions at short distances, but does still lead to the 
buildup of the Hall potential shown in Fig. 1 (c) In all 



our classical calculation, skipping orbits play only a mi- 
nor role for the transport and the frequent collisions with 
other electrons lead to a detachment of the trajectories 
from the edges. 

This observation poses the question whether in an 
interacting quantum-mechanical many-body calculation 
edge-states prevail. The electric Hall field and dissipa- 
tion at the contacts increase with increasing current and 
inelastic scattering events may provide another mecha- 
nism to change from an edge-transport picture to bulk 
transport diagonally across the device. The experimen- 
tal observations of the Hall potential similar to the one 



shown in Fig. 1(b) on a QHE plateau demonstrates that a 



quantized resistivity can occur hand-in-hand with the di- 
agonal transport picture. Transport along the edges is in 
our classical model not compatible with the cquipotential 
boundary conditions at the contacts, since the arrival of 
edge current in the contact would increase the electron 
density locally and lead to deviations of the potential 



from the prescribed value. If in the quantum-mechanical 
case a self-consistent calculation^ can remedy this situ- 
ation requires further investigations. 



VI. CONCLUSIONS AND OUTLOOK 

In conclusion, we have developed a highly efficient clas- 
sical particle-based simulation scheme for semiconduc- 
tor devices, which incorporates the spatial and temporal 
sequence of injection and removal events happening at 
Ohmic contacts. The classical Hall potential did emerge 
from our calculation as the time-averaged self-consistent 
potential. The choice of local boundary conditions within 
the metallic contacts has a crucial influence on the re- 
sulting global particle distribution inside a Hall bar. Our 
method extends previous models^ and can be readily 
adapted to include charged gates and different device ge- 
ometries. We expect that our new computational method 
opens a window towards a first-principles treatment of 
interaction effects in semiconductor devices and acceler- 
ates the development of the next generation of transport 
codes for realistic device settings. 
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